Magnetization profiles and NMR spectra of doped Haldane chains at finite 

temperatures 
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Open segments of S = 1 antiferromagnetic spin chains are studied at finite temperatures and fields 
using continuous time Quantum Monte Carlo techniques. By calculating the resulting magnetization 
profiles for a large range of chain lengths with fixed field and temperature we reconstruct the 
experimentally measured NMR spectrum of impurity doped Y2BaNii_ a; Mg :E 05. For temperatures 
above the gap the calculated NMR spectra are in excellent agreement with the experimental results, 
confirming the existence of S = 1/2 excitations at the end of open S — 1 chain segments. At 
temperatures below the gap, neglecting inter chain couplings, we still find well defined peaks in the 
calculated NMR spectra corresponding to the S = 1/2 chain end excitations. At low temperatures, 
inter chain couplings could become important, resulting in a more complicated phase. 

75.10.-b, 75.10 Jm, 75.40.Mg, 75.50.Ee 



I. INTRODUCTION 

The physics of low-dimensional spin systems is ex- 
tremely rich and shows many surprising effects. Using 
a valence bond picture it has been suggested that open 
segments of 5 = 1 spin chains should have S = 1/2 
excitations localized at the ends of the chain [1,2], that 
are experimentally observable. In a recent experiment 
by Tedoldi et al [3] the presence of these 5 — 1/2 chain 
end spins was beautifully demonstrated by measuring the 
NMR profile of Mg doped Y 2 BaNi0 5 , a model 5 = 1 spin 
chain system, showing well defined satellite peaks indica- 
tive of the chain end excitations. These satellite peaks 
were, perhaps surprisingly, very well defined at temper- 
atures close to and well above the gap, A ~ 100K. At 
temperatures below the gap, where one would expect the 
signal from the chain end spins to be most clearly observ- 
able, the satellite peaks were smeared into a single broad 
line. The inter chain coupling J- 1 in Y2BaNiOs is very 
small (| J x / J | < 5 x 1(T 4 , J ~ 285K) [3] and it therefore 
remains a question why no well defined satellite peaks 
were observed below the gap. Bulk susceptibility and 
magnetization measurements [4] on Zn doped Y2BaNiC>5 
have also been interpreted in terms of 5 = 1/2 chain end 
excitations, but in this case at temperatures all the way 
down to 4-10K, well below the gap. Very interestingly, 
a sub-Curie power law divergence was observed at tem- 
peratures below 4K hinting at the presence of a gapless 
phase. In light of these two experiments it is then natural 
to ask, from a theoretical point of view, if chain end ex- 
citations are observable at finite fields and temperatures 
well above the gap and when a simple one-dimensional 
model will fail to describe the experiments. In the present 
paper we systematically address these issues by calculat- 
ing the on-site magnetization, Sf, at finite temperatures 
and fields using continuous time quantum Monte Carlo 



techniques and reconstructing the associated NMR spec- 
tra for a given concentration of impurities. 

It is by now well established that integer spin chains 
display a gap above a spin liquid ground-state [5] as op- 
posed to half-integer spin chains which are gapless. The 
gap in integer spin chain systems has been observed ex- 
perimentally in numerous compounds [6] as well as in 
many numerical studies [7]. For S = 1 chains the gap is 
known to be A ~ 0.4105 J [8]. For integer spin chains, de- 
scribed by a simple antiferromagnetic Hcisenbcrg model, 
it is possible to understand the low temperature physics 
using the intuitive valence bond picture. From this point 
of view the ground-state is close to the valence bond solid 
(VBS) state. Each integer spin 5 is decomposed into 25 
spin 5 = 1/2 in the symmetric state and each of the 
5 = 1/2 on neighboring sites are paired into singlets. It 
is possible to determine Hamiltonians where this VBS 
state is the exact ground-state. For 5 = 1 spin chains 
the corresponding Hamiltonian is rather simple and apart 
from the usual Heisenberg term includes a biquadratic 
term -/3(Sj • S i+1 ) 2 . For (3 = -1/3, the AKLT point, 
the VBS state is the exact ground-state as was shown 
by Affleck-Kennedy-Lieb and Tasaki [9] (AKLT). Under 
periodic boundary conditions the VBS state is always a 
singlet. However, as was noted by Kennedy [1] for open 
5=1 chains, the ground-state is four- fold degenerate 
in the thermodynamic limit. In the VBS picture this is 
easy to understand since the 5=1/2 spins at the end of 
the chain are not paired into a singlet state, leading to 
a four-fold degeneracy for a 5 = 1 system. In general a 
spin-5 chain should then have spin-5/2 chain-end exci- 
tations. Many real compounds are well described by the 
antiferromagnetic Heisenberg chain: 



H — J Si ■ Si 



(1) 



Following the above discussion the ground-state should 
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then be close to the VBS state. For aS = l system, the 
low-energy states of an open segment is then a singlet and 
a triplet separated by a gap decaying exponentially with 
L, the length of the segment. For L odd the interaction 
between the two chain end excitations is ferromagnetic 
and the triplet is the lowest lying state, for L even the 
interaction is antifcrromagnetic and the singlet is lowest. 

In the following we will focus on S = 1 spin chains 
relevant to the experiments on Y 2 BaNi05. Numerically 
the 5 = 1/2 chain end excitations have been observed 
by various methods : exact diagonalisation [1] , Quantum 
Monte Carlo (QMC) [10], Density Matrix Renormaliza- 
tion Group (DMRG) [11]. They subsist all through the 
Haldane phase — 1 < /3 < 1 [11] but as (3 is increased 
from to 1 the peak in the onsite magnetization, 5f, 
moves away from the first site of the chain. At /? = 1 a 
dimerized phase is reached where the 5 = 1/2 chain end 
excitations cease to exist. At the AKLT point, (3 = —1/3, 
it can be analytically shown that the decay of the on-site 
magnetization away from the end of the chain has a pure 
exponential form with a length scale exactly equal to the 
bulk correlation length. Experimentally the S = 1/2 
chain end excitations have been observed in EPR mea- 
surements in S = 1 chains doped with non magnetic 
ions [2,12-15] and even in microscopically fractured un- 
doped systems [16]. Initially, experiments measuring the 
Schottky anomaly in Zn doped Y 2 BaNi05 [17] were in- 
terpreted in terms of S = 1 excitations but later work [18] 
showed that an explanation in terms of 5 = 1/2 chain 
end excitations is also possible when anisotropy is taken 
into account. These studies have mainly dealt with low 
temperature (i.e. T <C A) behavior of these systems 
where the only relevant states are the low-lying singlet 
and triplet. In the recent experiments by Tedoldi et al [3], 
a 89 Y NMR study of doped system YaBaNi^Mg^Os 
was performed, clearly indicating the 5=1/2 chain end 
excitations for temperatures well above the gap where 
many excited states apart from the low-lying triplet and 
singlet should be populated. In order to understand the 
NMR spectra in this high temperature range we simu- 
lated the isotropic Heiscnbcrg spin 1 chain using a con- 
tinuous time Quantum Monte Carlo algorithm at finite 
temperature and finite field. The QMC method sam- 
ples the partition function so that all (excited) states 
are correctly treated. For the purpose of understand- 
ing the physics at temperatures close to T ~ we have 
performed Density Matrix Rcnormalisation Group cal- 
culations of the NMR spectra at T = 0. Our results for 
temperatures above the gap are in very good agreement 
with the experimental data and clearly show that the 
5 = 1/2 chain end excitations indeed are observable at 
such high temperatures. At temperatures below the gap 
down to T = we find very well defined satellite peaks 
that at T = only are broadened by the distribution of 
chain lengths. In the experiments of Tedoldi [3] et al, 
no clear satellite peaks were observed at T — 77K, well 



below the gap. This could be due to a limited resolution 
at this temperature [19]. In the experiments of Payen et 
al [4] on doped Y 2 BaNi05 the results were interpreted 
in terms of the chain end excitations at temperatures 
down to 4K and it therefore remains unclear when and if 
higher dimensional physics become important in doped 
Y 2 BaNi0 5 . 

The paper is organized as follows. In section II we dis- 
cuss the numerical methods used. Section III presents 
our results for the on site magnetization for different tem- 
peratures and fields. In section IV we turn to a discus- 
sion of the NMR spectra obtained from the magnetiza- 
tion profiles and finally in section V we discuss possible 
explanations why signals from the 5=1/2 excitations 
are not clearly seen below the gap. 

II. NUMERICAL METHOD 

We use a single-cluster continuous time [20] version of 
the loop algorithm [21]. Spins 1 are simulated by two 
spins 1/2 with symmetrized boundary conditions [22,23] 
in the imaginary time direction. In order to take into ac- 
count the magnetic field, we use a hybrid method [24] 
where field-dependent loop flip probability and global 
flip of worldlines are used. Each Monte Carlo run con- 
sists of 10 7 sweeps (a measurement is taken every 10 
sweeps to avoid autocorrelations) preceded by 10 sweeps 
of thcrmalization. The DMRG [25] calculations were per- 
formed keeping m = 243 states for chain lengths between 
L=l,100 and m = 81 states for chain lengths L=101,501. 
Only negligible differences were observed between the re- 
sults for m = 81 and m = 243. The calculations were 
performed using the total z-component of the spin, S^ ot , 
and parity, P, with respect to a reflection about the mid- 
dle of the chain as quantum numbers. All the magneti- 
zation profiles were calculated in the 5 t z ot = 1,P = — 1 
state. 

III. MAGNETIZATION PROFILES 

A. Influence of temperature 

The first question we address is the temperature depen- 
dence of the on-site magnetization Sf . Due to symmetry 
the on-site magnetization is always zero in the absence of 
a magnetic field. In order to compare to experiments we 
therefore perform calculations in the sector with S^ ot = 1 . 
In Fig. 1 we show the magnetization profile for an open 
chain of 29 spins for various temperatures in the sec- 
tor 5 t z ot = 1. Filled circles are T=0 results obtained by 
DMRG, the other symbols are QMC data. We observe 
the well-known alternating structure due to the antifcr- 
romagnetic interactions which tends to be diminished as 
the temperature is raised. For the lowest temperature 
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figure 1: Magnetization profiles in the S^ ot — 1 for an 
open L=29 chain for different temperatures. The T=0 
profile (filled circles) is obtained by DMRG, other tem- 
peratures by QMC. 

T = 0.1J (well below the gap), the profile is almost iden- 
tical to the ground-state result. The decay of the on-site 
magnetization away from the end of the chain is at T = 
known to be exponential [26,11] with a length scale equal 
to the bulk correlation length £ ~ 6. At high temperature 
the staggered magnetization decays rapidly from its value 
at the end of the chain corresponding to a shorter corre- 
lation length. At the highest temperatures, the middle 
spins (around i = 15) have become polarized and show 
an average magnetization of < Sf >= \ ~ 0.034. At 
temperatures relevant to the experiments of Tcdoldi et 
al [3], T = 0.3 J - 0.8 J, the structure from the S = 1/2 
chain end excitations is still clearly visible. All these 
results confirm the existence of S = 1/2 excitations at 
temperature above the gap, even if the magnetization at 
the end of the chain is significantly reduced above the 
gap. 

B. Influence of magnetic field 

Lifting the constraint of working in the 5* t z ot = 1 sector 
we now consider the on-site magnetization as a function 
of the applied magnetic field at temperatures and fields 
relevant to the experiments. We calculated the depen- 
dence on the magnetic field of the profiles for two tem- 
peratures : T = 0.1J <C A and T = 0.526J > A. 

For T = 0.1J (figure 2), the zero-field magnetization 
is zero for all spins (by symmetry). At this temperature, 
the only magnetic states are in the sectors 5f ot = ±1 
of the triplet (which is the lower state for an odd chain) 
and are the only ones affected by the field. As we increase 



figure 2: Magnetization profiles for an open L=29 chain 
in different magnetic fields at T = 0.1J. 

the field, the Sf ot — 1 state will be energetically favored 
and hence the magnetization starts to take on the alter- 
nating structure showing the presence of the chain end 
excitations. At sufficiently large fields we should obtain 
results comparable to the results discussed in the previ- 
ous section obtained in the absence of a magnetic field 
but in the S^ ot = 1 sector. This is clearly the case and 
for H = 0.3J, we nearly recover the zero-field structure 
of the Sf ot = 1 sector at T = 0.1J shown in Fig. 1. If 
the field is too strong (say H > 0.3J), magnetic states 
from different sectors above the Haldane gap contribute 
to the thermal expectation of the on-site magnetization 
and, even if the alternating structure is still present, we 
don't recover the results from the zero-field case. This 
effect is visible on the H = 0.4 J data (left triangles). 
For this field, < > for odd spins (which all have a 
positive expectation value) is nearly equal to the value in 
the Sf Qt = 1 sector (filled circles) but for even spins the 
value for H = 0.4J is greater than in the S^ ot = 1 sector 
(this is particularly clear for the middle even spins). 

At higher temperature (T = 0.526 J), slightly above 
the gap A ~ 0.35 — 0.4 J (see figure 3), end effects are 
smaller, especially for low fields. At H = 0, the on-site 
magnetization is zero by symmetry and is not shown in 
the figure. Increasing the field, we notice two effects: the 
first is simply progressive magnetization of the middle 
spins, unaffected by end effects. When we later discuss 
the NMR spectra this will result in a positive shift of 
the NMR spectra for Y atoms close to the middle of the 
chain segments, i.e. of the central peak. The second 
effect is the progressive magnetization of the end spins 
with increasing field. The largest value of the field is 
equal to the field used in the experiments by Tcdoldi 
et al [3]. Clearly at this field and temperature there is 



3 



0.06 



0.04 



0.02 



0.00 



-0.02 



• H=0.08J 
i H=0.06J 

• H=0.04J 

• H=0.02J 



■ft 
if 



10 15 20 25 29 



figure 3: Magnetization profiles for an open L=29 chain 
in different magnetic fields at T = 0.526J. 

a signal from the S — 1/2 chain end spins. Since the 
temperature is large, the correlation length is small and 
the end effects are only visible for the first few spins away 
from the chain end. 



C. Estimation of correlation length 

As previously mentioned, the on-site magnetization 
should decay exponentially away from the end of the 
chain with a characteristic length related to the bulk 
correlation length. In order to compare to the experi- 
ments, we estimated the correlation length for different 
temperatures for the longest chains studied in a field of 
H = 0.08 J. We fit the exponential decay of the on-site 
magnetization away from the end of the chain using the 
following form: 



(Si) ((- 
+(s z h ). 



"UTJ 



+ (- 



L-i \ 



(2) 



Here, L is the length of the chain and the last term (Sf/) in 
this equation corresponds to the magnetization induced 
by the field at high temperature. Our results for the 
correlation length obtained using this method are shown 
in table I and are in agreement with previous numerical 
results [23] [27] and the experimental results taken from 
NMR experiments [3]. 



T(J) 


0.993 


0.702 


0.526 


0.27 




1.57 ± 0.23 


2.37 ± 0.37 


2.97 ± 0.13 


4.62 ± 0.11 



TABLE I. The temperature dependence of the correlation 
length as determined from the on-site magnetization profiles 
in a field of H = 0.08J. 



IV. NMR SPECTRA 

In the recent 89 Y NMR study by Tedoldi et al [3], 
of doped Y2BaNii_ a; Mg 2: 05 the Y is sitting roughly at 
an equal distance between two 5=1 chains. The non- 
magnetic Mg atoms will tend to break the chains into 
segments of different lengths and the NMR signal should 
then be proportional to the sum of the two nearest Ni 
spins sitting on different segments on different chains. 
Due to the rather large variation in chain lengths and 
the resulting variation in the on-site magnetization the 
satellite peaks should be broadened. Averaging over a 
large number of chain lengths it is then possible to re- 
construct the observed NMR spectra up to an overall 
rescaling factor. 

In order to perform this calculation wc have simulated 
open spin 1 segments in a magnetic field of H = 0.08J, 
corresponding to the experimental situation, for four 
temperatures between ~ 0.2J and ~ 1J equal to the 
experimental values, for chain lengths between L=l,200. 
At zero temperature we have computed the profiles in a 
slightly different fashion by selecting the S^ ot = 1 sub- 
space for chain lengths between L=l,501 using DMRG 
techniques. For very short even length chains, which are 
quite probable, this approach might lead to a slight bias 
since the magnetic field applied in the experimental situa- 
tion might not be sufficiently strong to select the 5* ot = 1 
state as the ground-state. However, we estimate this ef- 
fect to be rather small. 

Assuming that the distribution of the Mg impurities 
is uncorrelated we use a geometric distribution for chain 
lengths where the probability of having a chain of length 
L is given by: 



p(L) = x(l - x) 1 



(3) 



Here x is the impurity concentration. The mean chain 
length is then 



1=0 



1 

X 



1, 



(4) 



which, if x is sufficiently small, reduces to L = -. The 
non magnetic impurity concentration we used is x — 5%, 
which corresponds to a mean chain length of 19. One 
should note that the geometric distribution is quite wide 
and even though large chain lengths are less probable 
than short ones they contribute significantly to the ob- 
served NMR spectrum. The most probable chain length 
is in fact L = corresponding to an impurity site and 
the possibility of having one or more impurities next to 
each other is non-negligible. Hence, in order to obtain 
reliable results we have computed profiles for chains of 
long lengths. For T ^ 0, we calculated profiles for chains 
of lengths 1 up to 200, and for T = we simulated chains 
of up to 500 sites [28]. 
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According to the experimental situation, we randomly 
choose two lengths of chains and compute the sum of the 
two chains magnetizations site by site as we go along the 
chains. The two sites whose magnetizations are added 
at each step stand for the two next nearest neighbors 
spins of the NMR probe in the experiments [3]. When 
the end of one chain is reached we insert a site with zero 
on-site magnetization corresponding to the impurity. We 
then select a new chain following the geometrical distri- 
bution and continue to go along the system until the end 
of another chain is reached and we repeat the process 
till the thermodynamic limit is reached. If the length 
of the chain is we immediately insert an impurity site 
resulting in two or more impurities sitting next to each 
other. The NMR spectrum is then computed by plotting 
the histogram of the effective field seen on the Y site, i.e. 
the sum of the on-site magnetization coming from the 
two sites sitting closest to the Y on the two chains. This 
procedure should be quite close the experimental situa- 
tion. In order to compare to the experimental data we 
rescale the y and x axis of the experimental NMR spec- 
trum by a scale factor. This is done by adjusting to the 
first satellite peak (for the x axis we use a scale factor 
of roughly 0.8 — 0.9 which is in good agreement with the 
experimental estimate of hyperfine interaction) . 

In figures 4 and 5, we show the results of our calcula- 
tions for the five temperatures and compare them to the 
experimental spectra. The central peak corresponds to 
the situation where the effective field at the Y site can- 
cels either due to the fact that the two contributions are 
of roughly equal and opposite sign or because the Y is 
sitting close to the middle of both of the two chain seg- 
ments where the on-site magnetization is very close to 
zero. If the Y is closest to the first site of a chain seg- 
ment it should experience a positive shift away from the 
central peak and assuming that the magnetization on the 
other chain is close to zero we should observe the most 
pronounced satellite peak corresponding to the average 
value of the on-site magnetization on the first site. If we 
have exponentially localized excitations at the end of the 
chains the most probably value of the magnetization on 
the "other" chain is close to zero. However, this is of 
course not strictly the case, resulting in a broadening of 
the peaks. If the Y is close to the second site of one of 
the chain we should observe a satellite peak shifted in the 
negative direction away from the central peak but with 
a smaller absolute shift. It is then possible to identify 
the various satellite peaks with the on-site magnetiza- 
tion as one goes away from the end of the chain. As the 
distance from the end of the chain is increased the posi- 
tion of the satellite peak should move closer and closer 
to the central peak. Due to the significant broadening 
of the peaks only the first few satellite peaks are resolv- 
able. With exponentially localized excitations at the end 
of the chains there is also a small probability that the Y is 
close to firsts sites on both chains (for example first- first, 



first-third or second-second, second-fourth...). The- 
oretically, this should result in a number of secondary 
peaks beginning at twice the absolute value of the main 
peaks, also contributing to a general broadening of the 
spectrum. However, the intensity of such peaks should be 
rather small. In our calculations at T = 0, a secondary 
peak corresponding to the Y close to the firsts sites on 
both chains is clearly visible at twice the value of the first 
satellite peak (See Fig. 5). 
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figure 4: Calculated and experimental NMR spectra for 
an impurity concentration of x = 0.05 at T = 0.993J, 
0.702J, 0.526J and 0.270J. The results are obtained 
using QMC techniques at finite fields and temperatures 
for chain lengths of up to L = 200. Note that the last 
experimental spectrum has less resolution than at the 
other temperatures (the separation between the points 
corresponds to the experimental resolution) . 

The computed spectra are in extremely good agree- 
ment with the experimental results for the three higher 
temperatures. Note that the central peak is always 
greater in the NMR experiments than in our simulations, 
this might be due to a different distribution of impurities 
in the sample than the simple geometrical distribution 
(we don't consider correlations between impurities in this 
model) . There is a central peak for all the temperatures; 
the position of this peak is increased as temperature is 
raised due to the fact that (middle) spins are more easily 
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figure 5: Calculated NMR spectra for an impurity con- 
centration of x = 0.05 at T = 0. The results are obtained 
using DMRG techniques at zero temperature and zero 
field in the 5 t z ot = 1 subspace. Chain lengths of up to 
L = 501 are used in the calculation. The peak marked 1' 
corresponds to a site next to two chain ends correspond- 
ing to the secondary peaks discussed in the text. 

magnetized by the magnetic field for high temperatures, 
see section III B. The spectrum is wider at lower temper- 
atures since the magnetization profiles are much more de- 
veloped, as discussed in the previous section, resulting in 
larger shifts. For each figure, apart from the central one, 
several peaks can be observed. As explained, the shift 
from the central peak to the others correspond to the 
intensity of the magnetization on neighboring sites. On 
the figures, we have labelled the different satellite peaks 
according to the proximity to a site on the chain. As 
can be clearly seen, the neighboring magnetic structure 
of a site is more and more well defined as the tempera- 
ture is lowered. For example, at the highest temperature 
T = 0.993 J, only the two first satellite peaks are resolved, 
the more distant satellites cannot be distinguished and 
contribute to the central peak. For T = 0.526 J, up to six 
satellite peaks are resolved, and even more at the lowest 
temperature, T = 0.27J. This temperature is well below 
the gap and corresponds to the T = 77K experimental 
results in Ref . [3] . In this case the experiment show a sin- 
gle broad line and the signal from satellite peaks is lost. 
Tcdoldi et al [3] argues that one possible explanation of 
this is the large distribution of the chain lengths which 
could become more important at temperatures below the 
gap. This effect is taken correctly into account in our cal- 
culations and since we observe well-defined satellite peaks 
well below the gap and at T — the explanation is likely 
elsewhere. The resolution of the T = 77K spectrum is 
less than for the other temperatures [19] and it is possi- 



ble that higher resolution experiments will be capable of 
resolving the satellite peaks. 

Finally, in figure 5 we show results calculated at T = 
using DMRG techniques. The calculation is performed 
in zero field in the 5* ot = 1 subspace. Many well defined 
peaks are now visible. In particular we see secondary 
peaks corresponding to the case where the Y is between 
sites which both have large on-site magnetizations (!'). 



V. DISCUSSION 

Our results clearly indicate that the S = 1/2 chain 
excitations should be visible at temperatures well above 
the gap in a finite field. It is well known that the gap in 
S = 1 chains corresponds to a single magnon branch at 
k = 7T. Around k = a two-magnon continuum should 
begin at energies of 2A and a three-magnon continuum 
at k — 7r should begin at an energy of 3 A [29]. It is nat- 
ural to expect that chain end effects are visible until the 
two-magnon and three-magnon continuum become signif- 
icantly populated, i.e. for temperatures above 300K in 
agreement with the numerical and experimental results. 

The discrepancy between our theoretical results and 
the experimental results for temperatures below the gap 
could have several explanations. As suggested by Tedoldi 
et al [3] anisotropy effects could be important. However, 
the single ion anisotropy term D = 8K is rather small 
and it is not obvious that it would give rise to a large 
effect. Such a term would split the low- lying triplet into 
a doublet and a singlet and could likely cause a system- 
atic shift in the spectra rather than a complete smearing 
of the profile. Still another explanation could come from 
the small coupling between the chains. However, it would 
seem unlikely that this effect could be important at tem- 
peratures as high as T = 77K. Payen et al [4] argue that 
such effects become important at temperatures around 
4K, much lower than the gap, A ~ lOOif . However, 
these experiments are concerned with bulk susceptibil- 
ity and magnetization which might be a less sensitive 
probe. Another simple explanation is that below the gap 
so many secondary peaks appear that the experimental 
spectrum appears as a single broad line. This is not un- 
reasonable considering the results shown for T = 77K 
in figure 4. Lastly, the actual impurity concentration is 
not known to a very high precision and a slight modifica- 
tion in this concentration greatly modifies the observed 
spectra. High resolution NMR experiments at low tem- 
peratures will be needed to resolve these issues. 
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